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INFLUENCE OF SOIL AND TOPOGRAPHY ON THE 
COMPOSITION OF A TREE COMMUNITY IN A 
CENTRAL BRAZILIAN VALLEY FOREST 


J. R. R. Pinto!, A. T. OLIVEIRA-FILHO? & J. D. V. HAY? 


The purpose of the present study was to investigate the correlations between 
environmental variables related to the substrate (soil and topography) and the 
distribution of tree species in a valley forest in the Chapada dos Guimaraes National 
Park, Mato Grosso State, Brazil. The forest lies in the Cerrado Domain of Central 
Brazil where the dominant vegetation cover is cerrado (woody savanna). A survey of 
the tree community registered all live individuals with circumference at breast height 
(cbh) > 15cm found in eighteen 600m? plots (total area 1.08ha). The substrate variables 
used in the gradient analyses were obtained from a topographic survey and from 
analyses of the chemical and physical properties of soil samples. A principal 
components analysis of soil and topography variables and a canonical correspondence 
analysis of the species-environment relationships produced similar results, separating 
both the substrate variables and the tree species abundances, mainly according to the 
two types of bedrock, sandstone or slate, underlying their soils, and secondly to the 
three topographic sectors recognized: Streamside, Mid Slope and Upper Slope. The 
differences in soil fertility and texture (related to the bedrocks) and the soil water 
regime (related to both soil texture and topography) were probably the chief factors 
determining the distribution of tree species in the forest. 


Keywords. Brazil, Cerrado Domain, Chapada dos Guimarães, multivariate analysis, 
soil catena, tree species distribution, valley forest. 


INTRODUCTION 


Of the world's ecosystems, tropical forests are presently attracting much attention 
from ecologists, not only because of their species diversity but also because of their 
composition and complex structure. One of the components of this complexity is 
the heterogeneity of habitats, which influences the spatial distribution of plant spe- 
cies. Spatial heterogeneity in the physical environment is an important factor con- 
tributing to the commonly high tree species diversity of tropical forests, as variations 
in resource availability in both horizontal (particularly soil chemical and physical 
properties, and ground-water regime) and vertical (canopy layering, rooting zones) 
directions allow niche differentiation among tree species (Fowler, 1988; Terborgh, 
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1992). Development of the gap-phase dynamics theory for tropical forests (see 
Denslow, 1987 for a review) has made it clear that temporal heterogeneity of the 
environment (disturbance history) also plays a very important role in determining 
tree species distribution. Unfortunately, these findings have also played down, to a 
considerable extent, the notion that the substrate plays the major role in determining 
the horizontal patterns of tree species distribution (Clark et al., 1998). 

Tropical riverine forests are typical examples of high environmental-vegetational 
heterogeneity occurring over relatively small areas caused mainly by strong varia- 
tion in the substrate, particularly ground-water regime and soil chemical and physi- 
cal properties (Oliveira-Filho et al., 1994; Metzger et al., 1997; Hall & Harcombe, 
1998; van den Berg & Oliveira-Filho, 1999; Rodrigues & Nave, 2000; Rosales et al., 
2001; Botrel et al., 2002). Although cerrado (woody savanna) is the dominant vege- 
tation of the Cerrado Domain in Central Brazil, riverine forests are a very common 
element of the landscape because forests fringe nearly all water-bodies in the region. 
This forest network depends on high soil moisture throughout the year providing a 
suitable habitat for a large number of typical moist forest species despite the long 
dry season of the region (Oliveira-Filho & Ratter, 1995, 2002). Most of the striking 
variation of riverine forests in Central Brazil, in both physiognomy and floristic 
composition, results from variation in topography and drainage, together with soil 
properties (Oliveira-Filho et al., 1990, 2001; Felfili, 1995; Silva Júnior et al., 1996; 
Schiavini, 1997; Ribeiro & Walter, 1998; Rodrigues & Shepherd, 2000; Marimon 
et al., 2002). 

In some steep valleys of Central Brazil, riverine forests lie alongside mesophytic 
forests (instead of cerrado or grasslands) making up the so-called valley forests 
(Eiten, 1983, 1994; Oliveira-Filho & Ratter, 2002), many of which are sheltered in 
canyons and flanked by escarpments of sandstone plateaux. Since the descriptive 
accounts of Askew et al. (1971), Ratter et al. (1973, 1978) and Eiten (1975), Central 
Brazilian valley forests have received little attention, and, to our knowledge, there 
are no studies on their vegetation-environment relationships. In this paper we assess 
the correlations between environmental variables related to the substrate (soils and 
topography) and the distribution of tree species in a valley forest of the Chapada dos 
Guimarães. Our initial hypothesis was that variations in tree species abundances 
were chiefly related to elevation above the valley bottom. 


MATERIALS AND METHODS 
The study area 


The study area is a valley forest in the Chapada dos Guimarães National Park 
(PNCG), Mato Grosso State, Brazil, that lies at the bottom of the canyon below the 
Véu de Noiva (‘Bridal Veil’) waterfall, on the Coxipozinho stream at 15?2425"S 
and 55°50'25"W (Fig. 1). The scenic beauty of this area makes the Véu de Noiva 
waterfall the second most visited attraction of the Park (Macedo & Pinto, 1996). 
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Fic. 1. A, Position of the Chapada dos Guimaraes National Park (PNCG) in Mato Grosso 
State, Brazil. B, Surface grid of the canyon of the Véu de Noiva waterfall showing positions 
of the three sample transects (I, II, III). Lines on the surface grid are at 12.5m intervals. 


Altitudes range from 300-800m in the Park, and from 340-610m in the canyon. 
The regional climate is classified as Kóppens Aw and Cw types (IBAMA, 1994), 
both characterized by a dry winter, from May to September, and a rainy summer, 
from October to March. The mean annual precipitation ranges from 1800-2000mm 
(Oliveira-Filho & Martins, 1986). The Chapada dos Guimarães plateau is formed 
of Devonian sandstones and rests on the Pre-Cambrian slates that predominate in 
the surrounding lowlands, the Baixada Cuiabana (Oliveira-Filho & Martins, 1986, 
1991). The predominant soils in the region are dystrophic quartzitic sands at the 
top and the base of the plateau, dystrophic red-yellow latosols (oxysols) on the low- 
lands, and litholic soils formed over either sandstone or slate (MMA, 1982). Most 
soils in the forest studied are shallow; rocky outcrops are frequent and slopes 
are steep (20-50º); sandstone bedrock is replaced by slate at lower elevations. 
The predominant vegetation cover of the Park is cerrado, with its wide array of 
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physiognomic forms; seasonal marshes are found where the water table is periodic- 
ally near or at the surface; forests occur as narrow galleries, valley forests, or small 
patches of semideciduous forest (IBAMA, 1994). 


Survey of the tree community and environmental variables 


We surveyed the tree community in 18 permanent plots placed on catena transects 
aligned at right angles to the line of the watercourse and crossing both slopes to end 
at the foot of the escarpments on each side. We placed the transects in three sites 
along the canyon: Transect I at the inner end near the waterfall, Transect II in the 
middle and Transect HI near the mouth (Fig. 1). We marked out a total of 18 plots 
of 600m? (total area 1.08ha) at three topographic sectors in each of the six slopes 
crossed by the transects: Streamside, Mid Slope and Upper Slope. Plot dimensions 
were 10m x 60m in the Streamside sectors and 20m x 30m in the Mid Slope and 
Upper Slope; the shortest side was aligned parallel to the transect. We identified 
to species level all live trees in the plots with a circumference at breast height 
(cbh) > 15.0cm, measuring this with a tape and estimating total height with the help 
of a pole. 

We carried out a topographic survey along each of the three transects and the 
watercourse, and prepared a contour map and associated surface grid. Following the 
method used by van den Berg & Oliveira-Filho (1999), we calculated for each plot 
the topographic variable “vertical distance to the watercourse” from the difference in 
elevation between the plot centre and the valley bottom at the transect crossing. We 
collected four 0.5 litre soil samples from depths of 0-20cm at the midpoints between 
the corners and at the centre of each plot. The five soil samples were bulked and 
subsampled for each plot. Chemical and physical analyses were carried out at the 
Soil Laboratory of the Universidade Federal de Lavras following the procedures 
proposed by EMBRAPA (1997). 


Environmental analysis and tree community descriptive variables 


To identify the main variation patterns of the substrate in the area we performed 
a principal components analysis, PCA (Kent & Coker, 1992), of topography and 
soil variables using CANOCO 4 for Windows (ter Braak & Smilauer, 1998). The 
variables included: soil physical properties (proportions of sand, silt and clay), soil 
chemical properties (pH, extractable P, exchangeable bases [K*, Ca?*, Mg?* and AP*] 
and sum of bases [SB = K*+ Ca** + Mg?], proportion of soil organic matter, 
and vertical distance to the watercourse. Variables expressed as proportions were 
previously transformed by the expression arcsin (x^). 

The pattern arising from the PCA indicated two large groups of plots separated 
essentially by soil fertility and texture, reflecting whether the local bedrock was 
sandstone or slate. These two groups could be further split into five types of site 
based on the topographic sectors (Fig. 2), grouped as follows: 
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FIG. 2. Surface grids of sections of the canyon of the Véu de Noiva waterfall at the three 
sample transects showing the distribution of plots and the underlying bedrock. Lines on the 


surface grid are at 10m intervals. 


e Sandstone/Upper Slope (plots A, F, G, M, N, S) 


e Sandstone/Mid Slope (plots B, E, H, L) 
e Sandstone/Streamside (plots C, D) 

e Slate/Mid Slope (plots O, R) 

e Slate/Streamside (plots I, J, P, Q). 
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We compared soil variables among the five sites with analyses of variance, applying 
Tukey-Kramer tests to variables with significant differences (Zar, 1996). Variables 
expressed as proportions were previously transformed by the expression arcsin (x!?). 


Species—environment relationship 


We performed a canonical correspondence analysis, CCA (ter Braak, 1987, 1995), 
to investigate the relationships between species abundance and environmental 
(substrate) variables in the 18 plots using CANOCO 4 for Windows (ter Braak & 
Smilauer, 1998). The species abundance matrix consisted of the number of trees 
per plot, but included only the 39 species with 10 or more trees in the total sample. 
As recommended by ter Braak (1995), all abundance values were log-transformed 
before being analysed, as their distributions were skewed towards a few very large 
values. The matrix of environmental variables per plot initially included the same 12 
variables used in the PCA. The variable vertical distance to the watercourse was 
log-transformed for the same reason. In order to reduce the environmental variables 
to the most significant ones, we performed a preliminary CCA with forward selec- 
tion of environmental variables associated to Monte Carlo permutation tests for 
statistical significance. This procedure eliminated only the variable P. Additionally, 
we eliminated K*, Ca?* and Mg? because they were strongly interrelated (high 
redundancy indicated by variance inflation factor >20) and best summarized by 
the synthetic variable SB. The remaining eight variables used in the final CCA 
were vertical distance to the watercourse, pH, Al**, SB, organic matter, sand, silt 
and clay. The Monte Carlo permutation test (ter Braak, 1988) was performed to 
assess the significance of the correlations between the overall species abundance 
distributions and these eight environmental variables. 

To examine the effects of environmental variables in particular species popula- 
tions, we calculated Spearman’s rank correlation coefficients, and their significance 
(Zar, 1996), between species abundances in the 18 plots and two chosen substrate 
variables. The species were the same 39 used in CCA and the environmental vari- 
ables were SB and the vertical distance to the watercourse because they summarize 
soil chemical fertility and moisture, respectively. 


RESULTS 
Environmental analysis 


The results of the PCA are shown in the ordination diagram prepared with the first 
two components (Fig. 3). The contributions of the first four components to the total 
variation, expressed by the eigenvalues, are 0.616, 0.274, 0.074 and 0.029, respec- 
tively. The proportional accumulated variances explained by these components are 
61.9, 89.0, 96.4 and 99.3%. Therefore, the first two components summarize most 
variance contained in the environmental data. 
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Fic. 3. Ordination diagram produced by principal components analysis (PCA) of 11 envir- 
onmental variables (soil and topography) from the 18 sample plots in the valley forest of the 
Véu de Noiva waterfall. Each plot is identified by a letter (see Fig. 2) and given a particular 
symbol according to site as defined by underlying bedrock and topographic position. 
The arrows represent the main variation trend (direction and magnitude) of environmental 
variables. Distance = vertical distance to the watercourse. 


The first component separates the plots into two large groups, based essentially on 
soil fertility and texture, reflecting the two underlying bedrocks, sandstone and slate 
(Table 1). The Sandstone group is concentrated on the left side of the diagram and 
is strongly correlated with higher concentrations of Al** and proportions of sand. 
The Slate group, on the right side of the diagram, is strongly correlated with higher 
soil pH, higher concentrations of Ca**, Mg”, K* and SB, and higher proportions of 
silt and clay. The four plots of the Slate/Streamside sites are further from the centre 
of the diagram than those of the Mid Slope, suggesting internal variation among 
Slate plots based on fertility, distance to watercourse and organic matter. The 
second component separates the plots into groups corresponding to the topographic 
sectors and is strongly correlated with vertical distance to the watercourse and with 
soil organic matter. The two Sandstone/Streamside plots appear at the top of the 
diagram well away from the others probably because of their very low soil organic 
matter. 

The comparison of soil variables among the five sites discriminated by the PCA 
indicates that most significant differences occur between the soils derived from 
either sandstone or slate (Table 1). In general, sandstone-derived soils are more 
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TABLE 1. Environmental variables (1 topographic and 11 soil types) in the 18 sample plots. 
Figures are means + standard deviations of N plots or soil samples collected in the five types 
of site defined by underlying bedrock and topographic position (US= Upper Slope, 
MS=Mid Slope, SS=Streamside). Where F tests rejected the null hypothesis (P<0.05), 
means followed by different bold letters indicate significant differences in Tukey-Kramer tests 
(P « 0.05) 























Sand Silt Clay Organic matter 

Vertical —— 
Site N distance (m) (Yo) pH (H,0) 
Sandstone/US 6 38+163a 75+4.7a 17+3.9 8+21b 42+12a 4.5+0.32 b 
Sandstone/MS 4 22+13.8ab 72+5.7 ab 17+4.1 1142.5 ab 2540.5 ab 4.2+0.10b 
Sandstone/SS 2 4+40.7b 70+5.7 ab 18+5.0 12+0.7 ab 1.4+0.0b 4.5+0.35 b 
Slate/MS 2 25+28ab 60+8.5b 27 46.6 132.1 a 2.00.5 b 5.3 +0.00 a 
Slate/SS 4 83.2b 65+2.6b 23+2.6 12+0.0ab 2.6+0.5ab 5.5+0.19a 
ANOVAS F 5410 4.40 3.022 3.919 5.77 19.34 

P 0.011 0.018 0.054 0.027 0.067 <104 

Extractable Sum of 

P Kt Ca Mg* AP* bases 
Site N (mg dm-?) (cmol, dm) 
Sandstone/US 6 5.8+1.72 324+11.8¢ 0.7+0.46b 0.2+0.00¢ 3.0+0.98a 1.00.38 b 
Sandstone/MS 4 4.5+0.58 41418.7be 0.5+0.13b 0.24+0.00¢ 2.2+0.66ab 0.8+0.17b 
Sandstone/SS 2 4.0+1.41 3743.5¢ 0.5+0.14b 0.2+0.00¢ 1.8+0.2l ab 0.8+0.14b 
Slate/MS 2 6.00.00 101416.3a 0.9+0.07 ab 0.6+0.42b 0.940.64 be 1.7+0.42b 
Slate/SS 4 162+20.5 66+12.8ab 2.741.10a 1.0+0.10a 20.6-0.7c 3.9+41.10a 
ANOVAS F 147 11.42 2.22 29.41 8.53 17.30 

P 0.253 «10 «107 «10 0.013 «107^ 


acidic, contain higher amounts of exchangeable aluminium and lower amounts of 
exchangeable K*, Ca?* and Mg, and are coarser textured than slate-derived soils. 
SB is significantly higher only for Slate/Streamside soils. There is a significant 
increase in the levels of AP* and organic matter with increasing elevation within the 
Sandstone group. The levels of extractable P and silt did not differ significantly 
among sites. 


Tree community 


The comparison of physiognomic variables of the tree community among the sites 
(Table 2) indicates that tree density in the Sandstone sites was higher at the Stream- 
side and lower in the Mid and Upper Slope while basal area was higher in the 
Streamside and Mid Slope and lower in the Upper Slope. In the Slate sites, basal 
area was higher and density was lower in the Streamside than in the Mid Slope. 
Mean tree height within the Sandstone and Slate sites increased with decreasing 
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TABLE 2. Physiognomic variables for the total sample and the five types of site defined by 
underlying bedrock and topographic position 


No. of No. of Density Basal area dbh* Height* 

Site plots Area (ha) trees (trees ha-!) (m? ha^?) (cm) (m) 
Sandstone/Upper Slope 6 0.36 408 1133 20.825 118 82 
Sandstone/Mid Slope 4 0.24 262 1092 26.550 13.5 9.3 
Sandstone/Streamside 2 0.12 191 1592 25.400 11.6 9.8 
Slate/Mid Slope 2 0.12 206 1717 23.100 11.0 7.8 
Slate/Streamside 4 0.24 269 1121 24.829 12.7 84 
Total sample 18 1.08 1336 1237 23.749 12.2 8.5 
*Means. 


elevation and the mean diameter was highest at the Sandstone/Mid Slope and Slate/ 
Streamside sites. 

The 148 species registered in the sample plots with their quantitative parameters 
calculated for the total sample and the five sites defined by the PCA are given in 
Table 3. Among the 110 species sampled with more than one individual, Trichilia 
claussenii, Guarea kunthiana, Xylopia emarginata, Sebastiana membranifolia, 
Metrodorea stipularis, Mollinedia schottiana, Erythroxylum anguifugum and Piper 
amalago are exclusive to the Slate/Streamside plots. Tapirira guianensis, Vatairea 
macrocarpa, Qualea multiflora, Diptychandra aurantiaca, Myrciaria floribunda, 
Xylopia aromatica, Licania kunthiana, Vochysia haenkeana, Manihot tripartita, Ficus 
krukovii, Ocotea corymbosa and Cybianthus brasiliensis were found only in the 
Sandstone/Upper Slope plots. Myrsine lancifolia was found only in the Sandstone/ 
Mid Slope and Trichilia pallida in the Slate/Mid Slope. No species in that group was 
exclusive to the Sandstone/Streamside. The composition of the tree community is 
quite different among sites in terms of most abundant species. Of the 54 extracted 
from the top 15 of each site, only eight (14.8%) are common to both Sandstone and 
Slate sites. Among the three Sandstone sites, nine species are common to the Upper 
and Mid Slope, four to the Mid Slope and Streamside, and only two (3.7%) to the 
Upper Slope and Streamside. The two Slate sites share only four species. These 
results show that Slope and Streamside sites are quite distinct floristically for both 
the Sandstone and Slate areas. 


Species—environment relationships 


The results of the CCA are given in the ordination diagrams (Fig. 4) and the correl- 
ation matrix (Table 4). The eigenvalues of the first two CCA axes were 0.418 and 
0.228, respectively, indicating a moderately ‘long’ gradient for the first axis (= 0.5), 
i.e. there is considerable species turnover along the gradient (ter Braak, 1995), and 
a ‘short’ gradient for the second axis (<< 0.5), i.e. most species occur throughout 
the gradient and vary essentially in abundance. CCA axes accounted for a low 
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TABLE 3. Tree species recorded. Nt, number of trees; BA, total basal area; Rs, Spearman 
rank correlations with respective significance (P) between the species abundance (number of 
trees) and two environmental variables: sum of bases and vertical distance to the watercourse 
(N=18 plots). P-values<0.05 are underlined in bold. Species ranked by descending Nt. 
Correlations are given only for the 39 species with Nt 10. (The suffix ‘aceae’ is omitted from 


family names.) 


Species Families 

Siparuna guianensis Aublet Siparun. 

Ecclinusa ramiflora Mart. Sapot. 

Calyptranthes paniculata Myrt. 
Ruiz & Pav. 

Protium guianense subsp. Burser. 


pilosissimum (Engler) Daly 
Unonopsis lindmanii R.E.Fries Annon. 
Pseudolmedia laevigata Trécul Mor. 
Sloanea tuerckheimii Donn. Sm. Elaeocarp. 


Matayba guianensis Aublet Sapind. 
Trichilia claussenii CDC. Meli. 
Eugenia florida DC. Myrt. 
Inga vera Willd. Fab. 


Protium heptaphyllum (Aublet) Burser. 
Marchand 

Coussarea hydrangeifolia Rubi. 
(Benth.) Benth. & Hook.f. 


Miconia chartacea Triana Melastomat. 

Nectandra cuspidata Nees Laur. 

Ormosia arborea (Vell.) Harms Fab. 

Sacoglottis mattogrossensis Humiri. 
Malme 

Ocotea pomaderroides Laur. 


(Meisn.) Mez 


Guatteria conspicua R.E.Fries Annon. 
Miconia albicans Triana Melastomat. 
Miconia matthaei Naudin Melastomat. 
Bocageopsis mattogrossensis Annon. 


R.E.Fries 

Cheiloclinium cognatum (Miers.) Celastr. 
A.C.Smith 

Cecropia pachystachya Trécul Urtic. 


Ocotea aciphylla (Nees) Mez Laur. 

Tocoyena brasiliensis Mart. Rubi. 

Miconia punctata (Desr.) D.Don Melastomat. 

Maprounea guianensis Aublet Euphorbi. 

Chrysophyllum amazonicum Sapot. 
T.D.Penn. 

Tapirira guianensis Aublet Anacardi. 


Oenocarpus distichus Mart. Arec. 


Nt 


88 
65 
50 


13 


13 


BA 
(m?) 


0.310 
0.994 
0.486 


0.209 


0.534 
0.772 
1.300 
0.421 
0.533 
0.087 
0.242 
0.308 


0.171 


0.145 
0.510 
0.229 
0.335 


0.400 


0.276 
0.065 
0.081 
0.269 


0.132 


0.285 
1.081 
0.062 
0.158 
0.200 
0.229 


0.265 
0.493 


Sum of bases Vertical distance 


Rs 


0.02 
— 0.02 
—0.22 


—0.56 


0.44 
—0.01 
—0.35 

0.27 

0.63 

0.34 
—0.11 

0.39 


0.61 


—0.32 
—0.22 

0.13 
—0.37 


—0.29 


— 0.29 

0.09 
—0.50 
— 0.53 


P 


0.952 
0.950 
0.370 


Rs 


0.10 
— 0.23 
0.49 


—0.40 


—0.18 
— 0.46 
— 0.30 
— 0.40 
— 0.43 
— 0.61 

0.37 

0.02 


—0.21 


0.27 
0.15 
0.27 
0.29 


0.41 


—0.44 
0.33 
—0.52 
0.56 


—0.47 


— 0.09 
0.48 
— 0.23 
—0.47 
0.05 
0.20 


0.49 
0.44 


P 


0.682 
0.350 
0.038 


0.097 


0.486 
0.055 
0.222 
0.096 
0.074 
0.007 
0.127 
0.951 


0.398 


0.280 
0.543 
0.270 
0.241 


0.095 


0.069 
0.175 
0.028 
0.015 


0.048 


0.722 
0.043 
0.361 
0.047 
0.836 
0.420 


0.039 
0.066 
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TABLE 3. (Cont'd). 


Species 


Sclerolobium paniculatum Benth. 


Hymenaea courbaril L. 
Physocalymma scaberrimum 
Pohl 
Machaerium brasiliense Vogel 
Pouteria ramiflora Radlk. 
Jacaranda copaia (Aublet) 
D.Don. 
Sapium glandulosum 
(L.) Morong 
Ocotea elegans Mez 
Inga cayennensis Sagot 
Qualea multiflora Mart. 


Miconia longifolia (Aublet) DC. 


Licania apetala (E.Meyer) 
Fritsch 

Priogymnanthus hasslerianus 
(Chodat) P.S.Green 

Micropholis venulosa 
(Mart. & Eichler) Pierre 

Casearia arborea 
(L.C.Rich.) Urban 

Sebastiania membranifolia 
Müll.Arg. 

Matayba elaeagnoides Radlk. 

Guarea guidonia (L.) Sleumer 

Quiina rhitidopus Tul. 

Protium spruceanum 
(Benth.) Engler 

Alchornea glandulosa 
Poepp. & Endl. 

Connarus perrottetii 
(DC.) Planchon 

Virola sebifera Aublet 

Amaioua guianensis Aublet 

Miconia splendens 
(Swartz) Griseb. 

Elaeoluma glabrescens 
(Mart. & Eichler) Pierre 

Apuleia leiocarpa 
(Vogel) Macbr. 

Sloanea guianensis 
(Aublet) Benth. 

Bellucia grossularioides 
(L.) Triana 


Families 


Fab. 

Fab. 
Lythr. 
Fab. 
Sapot. 
Bignoni. 
Euphorbi. 
Laur. 
Fab. 
Vochysi. 
Melastomat. 
Chrysobalan. 
Ole. 
Sapot. 
Salic. 
Euphorbi. 
Sapind. 
Meli. 
Quiin. 
Burser. 
Euphorbi. 
Connar. 
Myristic. 


Rubi. 
Melastomat. 


Sapot. 
Fab. 
Elaeocarp. 


Melastomat. 


Sum of bases Vertical distance 


Rs P Rs P 
—0.09 0.711 0.32 0.195 
0.09 0.720 —0.30 0.220 
—0.13 0.611 0.20 0.436 
0.36 0.142 0.05 0.836 
—0.25 0.324 0.18 0.469 
— 0.60 0.009 —0.25 0.314 
0.50 0.036 —0.24 0.336 
— 0.29 0.236 —0.05 0.836 
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TABLE 3. (Cont'd). 


Species 


Tetragastris cerradicola Daly 
Inga marginata Willd. 
Guarea kunthiana A Juss. 
Chrysophyllum gonocarpum 
(Mart. & Eichler) Engler 
Cordia bicolor A.DC. 
Alibertia macrophylla K.Schum. 
Licania sclerophylla 
(Mart.) Fritsch 
Piper arboreum Aublet 
Allophylus edulis (A.St.-Hil.) 
Radlk. 
Urera baccifera (L.) Gaud. 
Licania hoehnei Pilger 
Bauhinia longifolia (Bongard) 
D.Dietrich 
Heisteria ovata Benth. 
Spondias mombin L. 
Copaifera langsdorffii Desf. 
Sparattosperma leucanthum 
(Vell.) K.Schum. 
Sorocea guilleminiana Gaudich. 
Xylopia aromatica (Lam.) Mart. 
Myrsine lancifolia Mart. 
Myrciaria floribunda 
(Willd.) O.Berg 
Hirtella glandulosa Sprengel 
Mabea fistulifera Mart. 
Cybianthus brasiliensis 
(Mez) Agostini 
Vismia guianensis (Aublet) Pers. 
Miconia nervosa (Smith) Triana 
Manihot tripartita Müll.Arg. 
Casearia gossypiosperma Briquet 
Buchenavia capitata 
(Vahl) Eichler 
Licania kunthiana Hook.f. 
Himatanthus lancifolius 
(Müll.Arg.) Woodson 
Tabebuia serratifolia 
(Vahl) Nichols 
Schefflera morototoni 
(Aublet) Maguire, 
Steyerm. & Frodin 
Metrodorea stipularis Mart. 


Families 


Burser. 
Fab. 
Meli. 
Sapot. 


Boragin. 
Rubi. 


Chrysobalan. 


Piper. 
Sapind. 


Urtic. 


Chrysobalan. 


Fab. 


Olac. 
Anacardi. 
Fab. 
Bignoni. 


Mor. 
Annon. 
Myrsin. 
Myrt. 


Chrysobalan. 


Euphorbi. 
Myrsin. 


Clusi. 
Melastomat. 
Euphorbi. 
Salic. 
Combret. 


Chrysobalan. 


Apocyn. 
Bignoni. 


Arali. 


Rut. 


Nt 


ON O00 0 


BA 
(mr) 


0.099 
0.102 
0.139 
0.161 


0.372 
0.018 
0.023 


0.023 
0.025 


0.046 
0.055 
0.072 


0.120 
0.215 
0.301 
0.813 


0.013 
0.027 
0.028 
0.029 


0.085 
0.108 
0.006 


0.008 
0.009 
0.010 
0.014 
0.020 


0.022 
0.039 


0.041 


0.094 


0.106 


Sum of bases Vertical distance 
Rs P Rs P 
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TABLE 3. (Cont'd). 


BA Sum of bases Vertical distance 











Species Families Nt (ny) Rs P Rs P 
Diptychandra aurantiaca Tul. Fab. 3 0.132 
Xylopia emarginata Mart. Annon. 3 0.241 
Vatairea macrocarpa Fab. 3 0.260 


(Benth.) Ducke 
Erythroxylum anguifugum Mart. Erythroxyl. 2 0.004 




































































Ocotea corymbosa (Meisn.) Mez Laur. 2 0.006 

Mollinedia schottiana Monimi. 2 0.007 
(Sprengel) Perkins 

Piper amalago L. Piper. 2 0.007 

Ficus krukovii Standley Mor. 2 0.009 

Trichilia pallida Swartz Meli. 2 0.010 

Duguetia echinophora R.E.Fries | Annon. 2 0.020 

Vochysia haenkeana Mart. Vochysi. 2 0.032 

Gomidesia hebepetala Myrt. 2 0.038 
(DC.) O.Berg 

Calyptranthes strigipes O.Berg = Myrt. 2 0.039 

Aspidosperma spruceanum Benth. Apocyn. 2 0.080 

Sterculia apetala (Jacquin) Malv. 2 0.121 
H.Karst. 

Pera glabrata (Schott) Poepp. ^ Euphorbi. 2 0.127 

Attalea phalerata Mart. Arec. 2 0.186 

Aegiphila triantha Schauer Lami. 1 0.002 

Cybianthus guyanensis Myrsin. 1 0.002 
(A.DC.) Mig. 

Eugenia uniflora L. Myrt. 1 0.002 

Heteropterys byrsonimifolia Malpighi. 1 0.002 
A.Juss. 

Inga heterophylla Willd. Fab. 1 0.002 

Mouriri acutiflora Naudin Memecyl. 1 0.002 

Rhamnidium elaeocarpum Rhamn. 1 0.002 
Reissek 

Schefflera vinosa Arali. 1 0.002 


(Cham. & Schltdl.) 
Frodin & Fiaschi 


























Banara tomentosa Clos Salic. 1 0.003 

Endlicheria paniculata Laur. 1 0.003 
(Sprengel) Macbr. 

Heisteria densifrons Engler Olac. 1 0.003 

Roupala brasiliensis Klotzsch Prote. 1 0.003 

Aegiphila lhotskiana Cham. Lami. 1 0.004 

Apeiba tibourbou Aublet Maly. 1 0.004 

Ceiba speciosa (A.St.-Hil.) Maly. 1 0.004 
Ravenna 

Miconia minutiflora Melastomat. 1 0.004 


(Bonpl.) DC. 
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TABLE 3. (Cont'd). 


BA Sum of bases Vertical distance 




































































Species Families Nt (mn?) Rs P Rs P 
Ziziphus mistol Griseb. Rhamn. 1 0.004 
Chloroleucon tortum Fab. 1 0.005 
(Mart.) Pittier 
Virola albidiflora Ducke Myristic. 1 0.005 
Astrocaryum aculeatum Meyer — Arec. 1 0.006 
Eupatorium laeve DC. Aster. 1 0.006 
Ficus insipida Willd. Mor. 1 0.006 
Hyeronima alchorneoides Phyllanth. 1 0.006 
Fr.Allem. 
Acacia polyphylla DC. Fab. 1 0.007 
Ficus maxima P.Miller Mor. 1 0.007 
Cybianthus cuneifolius Mart. Myrsin. 1 0.008 
Terminalia glabrescens Mart. Combret. 1 0.009 
Astronium fraxinifolium Schott  Anacardi. 1 0.011 
Byrsonima coriacea (Swartz) DC. Malpighi. 1 0.015 
Maclura tinctoria (L.) D.Don. Mor. 1 0.019 
Styrax pohlii A.DC. Styrac. 1 0.026 
Antonia ovata Pohl Logani. 1 0.027 
Myrcia tomentosa (Aublet) DC.  Myrt. 1 0.048 
Guazuma ulmifolia Lam. Maly. 1 0.078 
Emmotum nitens (Benth.) Miers Icacin. 1 0.088 
Poecilanthe parviflora Benth. Fab. 1 0.115 
Simarouba amara Aublet Simaroub. 1 0.126 
Eriotheca gracilipes Maly. 1 0.255 
(K.Schum.) A.Robyns 
Totals: 148 52 1336 25.66 


cumulative percentage variance of species abundance data, 20.1% and 31.0% for 
axes 1 and 2, respectively, indicating that considerable “noise” remained unexplained. 
However, low percentage variances for species data are normal in vegetation data 
and do not impair the significance of species-environment relations (ter Braak, 
1988). In fact, CCA produced quite high values for both species-environment 
correlations (axis 1 20.97, axis 2=0.93). In addition, the Monte Carlo permutation 
test demonstrated that the species abundances and environmental variables were 
significantly correlated (F-ratio = 2.74, P « 0.01). 

The two underlying bedrocks appear as very distinct clusters of plots in the CCA 
diagram (Fig. 4A). The first canonical axis is most strongly (and negatively) correl- 
ated with Al, followed by sum of bases and pH (both positively), sand (negatively), 
and silt, clay and organic matter (all positively) (Table 4). These soil variables are 
also highly interrelated (see weighted correlation matrix in Table 4). The second 
canonical axis is more strongly influenced by the vertical distance to the watercourse 
while the other seven variables yielded very low correlations. The elevational 
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Fic. 4. Ordination diagrams produced by canonical correspondence analysis (CCA) of the 
number of individuals of the 39 most abundant species (V>10) sampled in eighteen 600m"? 
plots in the valley forest of the Véu de Noiva waterfall, and their correlation with environ- 
mental variables. Ordination of plots and environmental variables (A) and ordination of 
species (B) on the first two CCA axes are given in separate diagrams for clarity. Each plot is 
identified by a letter (see Fig. 2) and given a particular symbol according to site as defined by 
underlying bedrock and topographic position. The arrows represent the main variation trend 
(direction and magnitude) of environmental variables. Distance = vertical distance to the 
watercourse. Full species names are given in Table 3. 
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TABLE 4. Canonical correspondence analysis (CCA): intraset correlations in the first two 
ordination axes, and matrix of weighted correlations for the eight environmental variables 
supplied. Distance = vertical distance to the watercourse. Correlations with absolute values 
>0.5 are underlined in bold 


Intraset correlations Environmental variables 
Environmental Sum of Organic 
variables Axis 1 Axis 2 Distance pH Al bases matter Sand Silt Clay 
Distance — 0.48 0.76 - 
pH 0.75 0.32 0.14 — 
Al — 0.89 0.01 0.50 — 0.56 — 
Sum of bases 0.77 — 0.01 0.26 0.80 —0.49 = 
Organic matter —0.56 0.26 0.66 — 0.07 0.75 0.03 — 
Sand — 0.66 —0.01 0.29 —0.50 0.39 —0.45 0.26 - 
Silt 0.58 0.08 0.08 0.53 —0.26 0.50 —0.00 —0.95 - 


Clay 0.57 —0.14 0.62 0.22 —0.52 0.15 —0.68 —0.71 0.45 - 


gradient appears as a diagonal in the CCA diagram (Fig. 4A). These trends are also 
highly correlated with the differences in soil properties among sites given in Table 1. 

The species ordination by CCA (Fig. 4B) indicates correlation with the environ- 
mental variables and therefore suggests their preference for the five habitats defined 
by soil fertility and texture and topography. For instance, Trichilia claussenii, 
Cheiloclinium cognatum, Eugenia florida, Matayba guianensis and Coussarea 
hydrangeifolia show a preference for the Slate/Streamside site, while Protium 
heptaphyllum, Cecropia pachystachya, Sapium glandulosum, Machaerium brasiliense, 
Tocoyena brasiliensis, Unonopsis lindmanii, Hymenaea courbaril and Ecclinusa 
ramiflora prefer the Slate/Mid Slope sites. Physocalymma scaberrimum, Chryso- 
phyllum amazonicum, Inga vera, Ocotea aciphylla, Calyptranthes paniculata, 
Sclerolobium paniculatum, Siparuna guianensis, Pouteria ramiflora, Nectandra 
cuspidata, Maprounea guianensis, Sloanea tuerckheimii, Ocotea elegans, Protium 
guianense subsp. pilosissimum and Jacaranda copaia show a preference for the 
Sandstone/Mid Slope sites. Miconia punctata, Pseudolmedia laevigata, Guatteria 
conspicua and Miconia matthaei seem to prefer the Sandstone/Streamside, and 
Tapirira guianensis, Miconia albicans, Oenocarpus distichus, Ormosia arborea, 
Sacoglottis mattogrossensis, Bocageopsis mattogrossensis, Ocotea pomaderroides and 
Miconia chartacea, the Sandstone/Upper Slope sites. 

The species ordination by CCA (Fig. 4B) strongly agrees with their correlation 
coefficients with the two chosen substrate variables (Table 3), vertical distance to 
the watercourse and total bases. Of the 39 species involved in the analysis, 12 
showed significant correlations with either of the two variables and three with 
both variables. The three species with positively significant correlations with total 
bases, Trichilia claussenii, Coussarea hydrangeifolia and Sapium glandulosum, also 
appeared on the left side of the ordination diagram. In contrast, the six species 
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with negative correlations with total bases, Protium guianense subsp. pilosissimum, 
Miconia matthaei, Bocageopsis mattogrossensis, Ocotea aciphylla, Maprounea 
guianenis and Jacaranda copaia, appeared on the right side of the diagram. The four 
species with positively significant correlations with vertical distance to the water- 
course, Calyptranthes paniculata, Bocageopsis mattogrossensis, Ocotea aciphylla 
and Tapirira guianensis, are all concentrated in the upper part of the diagram, 
whilst the five negatively correlated with vertical distance, Pseudolmedia laevigata, 
Eugenia florida, Miconia matthaei, M. punctata and Cheiloclinium cognatum, are 
concentrated at the bottom. 


DISCUSSION 


Catenary variations of soil properties in tropical forests have commonly been attrib- 
uted to topography itself, because variations in steepness and different topographi- 
cal levels on a slope may determine, to a great extent, different ground-water 
regimes and these, in turn, usually affect soil texture and nutritional status 
(Bourgeron, 1983; Newbery & Proctor, 1984; Furley, 1992; Silva-Júnior, 1998; Silva 
Junior et al., 1996; Oliveira-Filho et al., 1997, 2001; van den Berg & Oliveira-Filho, 
1999). This was our initial model for studying the tree community gradients in 
this area. However, an unexpected pattern related to soil origin was revealed in 
the analyses and proved to be the leading factor influencing the distribution of tree 
species in the forest. A gradient related to topographic levels was also shown but 
was less pronounced. Both the environmental analysis (PCA) and the species— 
environment analysis (CCA) produced similar results and indicated that five sites 
could be defined for the valley forest based on underlying bedrock and topographic 
site. 

The influence of the local bedrocks in pedogenesis certainly caused the strong 
variation in soil texture and nutritional status within the canyon. Siltstones, slates, 
and other poor, pellitic rocks in Brazil are normally poor in Ca but rich in K and 
give rise to fine-textured soils; sandstones, on the other hand, often produce very 
dystrophic and coarse-textured soils (Resende et al., 1995). Previous studies have 
demonstrated that the strong dichotomy between soils originating from either slate 
or sandstone in the region of the Chapada dos Guimarães plateau and Baixada 
Cuiabana lowlands is reflected in the composition of the woody plant communities 
of both arboreal savannas ('cerradão”) (Oliveira-Filho & Martins, 1986, 1991) and 
gallery forests (Oliveira-Filho et al., 1990). In both cases, the authors identified soil 
chemical fertility and water storage capacity as the main factors determining differ- 
ences in species composition. In the case of arboreal savannas, soil fertility is the 
chief factor resulting in two well-known physiognomies found throughout Central 
Brazil: the ‘mesotrophic cerradão” and the ‘dystrophic cerradão” (Ratter, 1971, 1992; 
Ratter et al., 1973, 1978, 2003). Soil fertility was found to be the leading factor 
underlying the differentiation between adjacent evergreen forest on dystrophic soils 
and semideciduous forest on mesotrophic soils (Ratter et al., 1973). 
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In addition to differences in soil fertility, the two bedrocks give rise to soils with 
rather different water regimes. Sandstone-derived soils are much coarser and, there- 
fore, have stronger drainage and lower water storage capacity than slate-derived 
soils (Resende et al., 1988). It is quite difficult to separate the effects of ground- 
water and soil fertility on the differentiation of the plant communities growing on 
either soil. The importance of ground-water regime in species distribution in the 
valley forest is reinforced by the second level gradient dictated by the topographic 
levels. A single soil water gradient for the valley forest sites can be deduced, from 
highest to lowest water availability: 


Slate/Streamside — Sandstone/Streamside — Slate/Mid Slope — Sandstone/Mid 
Slope — Sandstone/Upper Slope 


Habitat differentiation also influenced tree community physiognomy. The concen- 
tration of large trees in the Sandstone/Mid Slope, Sandstone/Streamside and Slate/ 
Streamside sites may be due to stronger competition for light in the lower parts of 
the canyon and to higher water availability (Oliveira-Filho et al., 2001). Mean tree 
stature is less in upper sites, perhaps due to lower nutrient and water availability 
factors that may restrict tree growth (Lleras & Kirkbride, 1978). The higher density 
of smaller trees in the Slate/Mid Slope sites is probably related to soil depth, since 
the slate bedrock outcrops or is near the surface over the whole site. In these condi- 
tions, root systems are shallow, resulting in poor stability of the larger trees, which 
frequently fall forming gaps that provide conditions for the establishment of many 
small individuals. 

The present study is of an exploratory nature and the indication of habitat prefer- 
ences by particular species can only be confirmed after repeated studies or experi- 
mental work. For some species, there is already published evidence that the patterns 
shown are repeated elsewhere. Examples are the association of Pseudolmedia 
laevigata and Cheiloclinium cognatum with moister soils (Oliveira-Filho & Ratter, 
1995) and Trichilia claussenii with richer soils (Oliveira-Filho et al., 1997). For 
some species however there is always the possibility that environmental variables not 
directly related to the substrate may be the actual leading factors. Light regime is 
the major candidate through both gap-phase dynamics and light gradients caused 
by forest physiognomic variations (Denslow, 1987). The higher abundance of the 
pioneer tree Cecropia pachystachya in the Slate/Mid Slope, for example, probably 
results from the higher gap formation rate in those sites, as discussed above, and not 
from local soil fertility and moisture. In addition, canopy openness visibly increased 
away from the valley bottom and this certainly affects tree establishment in the 
understorey. Because light and topography gradients are related in this case, it is 
very difficult to differentiate their effects on plant distribution. The higher abun- 
dance of Miconia albicans in the Sandstone/Upper Slope sites, for example, may be 
due both to its preference for strongly drained sites and to its light-demanding 
behaviour (Espirito-Santo et al., 2002). To conclude, there is strong evidence that 
variations in the substrate, particularly soil moisture and fertility, were indeed 
important factors influencing tree species distribution in the valley forest. On the 
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other hand, there was evidence that other factors, particularly light, may also be 
involved. Great care should be taken when interpreting apparent substrate-species 
relationships. 
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